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Abstract 



. Sound scattering by a finite width beam on a single rigid body rotation 



vortex flow is detected by a linear array of transducers (both smaller than a 
flow cell) , and analyzed using a revised scattering theory. Both the phase and 
amplitude of the scattered signal are obtained on 64 elements of the detector 



> 

. array and used for the analysis of velocity and vorticity fields. Due to aver- 

o 

^ aging on many pulses the signal-to-noise ratio of the phases difference m the 

^ . scattered sound signal can be amplified drastically, and the resolution of the 



method in the detection of circulation, vortex radius, vorticity, and vortex lo- 
cation becomes comparable with that obtained earlier by time-reversal mirror 



^ ! (TRM) method (P. Roux, J. de Rosny, M. Tauter, and M. Fink, Phys. Rev. 

H : 

Lett. 79, 3170 (1997)). The revised scattering theory includes two crucial 
steps, which allow overcoming limitations of the existing theories. First, the 
Huygens construction of a far field scattering signal is carried out from a signal 
obtained at any intermediate plane. Second, a beam function that describes a 
finite width beam is introduced, which allows using a theory developed for an 
infinite width beam for the relation between a scattering amplitude and the 
vorticity structure function. Structure functions of the velocity and vortic- 
ity fields deduced from the sound scattering signal are compared with those 



obtained from simultaneous particle image velocimetry (PIV) measurements. 

Good quantitative agreement is found. 
PACS numbers: 43.30+m,43.35+d,47.32-y 
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I. INTRODUCTION 



Ultrasound wave scattering provides a non-perturbative tool to probing a hydrodynamic 
velocity field. Propagating sound waves are modified by their interaction with the flow that 
results in their amplitude and phase distortions and scattering at different angles [1]. The 
sound scattering by a flow was studied for the last 50 years theoretically [2-5] as well as 
experimentally [6-8] with a goal to develop a reliable technique to directly study dynamics 
of vorticity and velocity fields. Being developed and being reliable this tool can provide 
rather unique possibility to get information about dynamics of turbulence. However, in 
spite of the efforts the goal to measure the dynamical structure factor of the vorticity in a 
turbulent flow remains rather far from reach. In this paper we attempt to address much less 
ambitious but rather important experimental problem: is it visible to obtain reliable infor- 
mation about velocity and vorticity field of an axisymmetric large vortex either stationary 
or time-dependent from a resulting scattered signal of a finite width sound beam by a finite 
width receiver and rather close to the scattering region? 

The full structure of the scattering waves (amplitude and phase) can be related to a struc- 
ture function of a vorticity via analytical relation only in the Born approximation and in a 
far- field limit [4,5]. This theory however considers a beam, a sound emitter, and a receiver 
of an infinite width. Only Ref. [9] discusses an influence of a finite width of a Gaussian 
sound beam on scattering from a point vortex in a far- field limit. Finite width of these con- 
stituents particularly if they are smaller than a flow region, leads to additional diffraction 
effects, which are mixed with interference patterns resulted from the scattering and refrac- 
tion. This unavoidable nuisance effects should be unravelled from the intrinsic scattering 
signal, from which the vorticity could be determined. The aim of the studies presented in 
the paper is to develop a kind of analysis, which can allow determining the vorticity and 
velocity field without limitations described and to apply it to the experimental data on the 
sound scattering from a single either stationary or time-dependent vortex. 
The majority of the experiments in the past [10-12] were concerned with sound propagation 
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in turbulent flows and determination of the scattering cross-section, which can be compared 
with theoretically predicted one. Such approach cannot provide unambiguous answer about 
the scattering by a turbulent velocity field. Later on it was realized that the scattering 
waves may be useful for remote probing of turbulent flows [4,5]. Using this idea the scatter- 
ing problem was reformulated as the inverse problem to determine an instantaneous velocity 
field within a scattering region (laminar or turbulent) from measurements of the acoustic 
scattering at different angles similar to the light scattering problem. It was shown that in the 
approximation of a plane sound wave propagating through a velocity field with a continuous 
distribution of vorticity in a finite scattering domain there is a linear relation between the 
Fourier component of the scattering sound wave amplitude and the spatial and temporal 
Fourier transform of the vorticity component normal to the plane of the wave propagation 
in a far-field limit [4,5]. 

During the years, several experiments were conducted to study either stationary or time 
dependent vortex flows by sound scattering. Earlier experiments by Gromov et al [6] and 
later and more precise by Baudet et al [7] investigated the von Karman vortex street behind 
a cylinder at low Reynolds numbers. Re. The experiments showed rather good qualitative 
agreement with the theoretical predictions. Later on the experiment by Oljaca et al [8] 
on the ultrasound scattering from a stationary random vortex flow formed along the axis 
of a swirling jet at Re — 10^, showed a good quantitative agreement between measured 
and computed from separate velocity fleld measurements of the scattering sound amplitude. 
The main reason for the success was a satisfaction of rather strict limitations of the the- 
ory, which leads to the linear analytical relation between Fourier transform of the vorticity 
and the scattering acoustic signal [4,5]. In particular, it is required for an emitter to be 
perfectly flat and wider than the flow region, and for a detector to be positioned in a far 
fleld region. In the experiment [8] the characteristic size of the jet was about the sound 
wavelength and much smaller than the receiver and the sound beam widths. Hence, the 
requirements were perfectly met for the applicability of the theory. In this connection one 
can mention also Refs [13,14] where just structure functions of vorticity in a frequency do- 
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main were measured. Another approach was undertaken in a series of papers by Fink and 
collaborators [15-18], where the velocity field and vortex dynamics were reconstructed from 
a phase shift measured by a transducer array. Time-reversal mirror (TRM) method is used 
to amplify the effect of the flow on the sound phase distortion leading to a high resolution 
in the velocity detection. The acoustic technique is based on a geometrical approximation, 
which is justified for the fiow field and sound characteristics used in the experiments. In 
a forward scattering geometry the phase shift possesses the same information as a sound 
scattering but only in the strictly geometric acoustics limit. Diffraction processes, which 
are relevant for a finite wave length as well as finite width receiver and sound beam, require 
deeper understanding of sound scattering particularly in a projection of an acoustic signal 
from a near-field to a far-field region. On the other hand, there was a very recent attempt 
to achieve a goal similar to the formulated above, by using direct numerical simulations to 
analyze the whole structure of the scattered wave (amplitude and phase) and to compare it 
with the experimental measurements of the sound scattering on a vortex in a turbulent fiow 
[19]. The agreement found is not great, particularly for the amplitudes. Our approach is 
rather different: it is straightforward in acquisition of the phase and amplitude and simple 
to apply especially at higher frequencies. 

II. EXPERIMENTAL SET-UP AND MEASURING TECHNIQUES 

The experiments were performed in a fiow, which was produced inside a closed cylinder 
either with one rotating disk or in the gap between two coaxial either co- or contra-rotating 
disks (the so-called von Karman swirling fiow). As detection techniques we use ultrasound 
scattering and particle image velocimetry (PIV) to directly measure velocity field. 
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A. Hydrodynamic set-up 



The hydrodynamic cell consists of a vertical Plexiglass (Perspex) cyhnder of 29 cm inside 
diameter and 32 cm in height (see Fig.l). A swirling water flow is produced between two 
rotating plates driven independently by two motors with a maximum continuous torque of 
13Nm. The motors are brushless sinusoidal ones controlled by velocity mode drivers via 
optical encoders. Each shaft is centered and is fixed by two bearings and sealed with Uniten 
rotating rings. The base of the device is suspended on shock absorbers and the shafts are 
isolated from vibration by the motor couphngs (Rotex GS designed for servo systems). The 
acoustic components (emitters and a receiver array) are mounted flush to avoid obstacles 
to the flow (see Fig. 2). The cell is completely flUed with deionized water with addition of 
lOppm of surfactant. Small amount of surfactant facilitates removing of air bubbles from 
water. Several steering conflgurations have been used: (i) Two steering plates of 100 mm 
diameter are positioned 100 mm apart, each plate has four blades of 20 mm width and 5 
mm thickness; (ii) the same kind of plates with 40mm diameter and 4mm blades thickness; 
(iii) Upper plate is a metal disk of 280 mm diameter, a lower plate is absent, and only 25 
mm diameter shaft is present with upper plane located at 50 mm from the bottom of the 
cell; (iv) Two elongated shafts of 25 mm diameter are set 86 mm apart, grinded at the edge 
to a shape of one blade of 12 mm width and 6 mm thickness. 

B. Principle of measurements 

1. Ultrasound transducers 

The measurement plane where the emitter and the detector array are placed, is chosen 
at the middle plane between the disks. The sound detects the velocity field in the plane, or 
complimentarily, the vorticity field perpendicular to that plane. Ultrasound pulses are sent 
by a specially designed emitter (see below) and detected by a hnear array of 64 acoustic de- 
tectors with 1mm spacing and 64xl0mm overall active area (from Blatek). The acquisition 
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system is built in a heterodyne scheme meaning that 64 lock-in amplifiers are utilized on the 
incoming analog signals [20]. The outgoing low frequency signals are passed to sample and 
hold components, which integrate the signal of each pulse and hold the value to be scanned 
by two PC acquisition cards (National Instruments pci6071) recording 128 channels, hence 
providing amplitude and phase components of the signal for 64 acoustic channels. The lock- 
in amplifiers were compared to a SR844 RF lock-in amplifier (Stanford Research Systems) 
and were found linear and reliable. The acquisition cards also function as controllers of 
timing with precision of 0.05 /is. The detector and electronics are optimized for 5.5MHz. 
Each of 64 preamplifiers is based on VCA2612 chips from Texas Instruments, and it is built 
directly around the detector. The entire circuitry consists of about 8,000 components. A 
special care was taken to block electromagnetic interference emitted from the motors: (i) 
installing filters on the motors and the power supply; (ii) the circuitry and the emitter cable 
are shielded in a triplex scheme where the outer shell is connected to chassis ground, the 
intermediate shell connected to a one point analog ground. 

We are able to measure amplitude and phase of the sound pressure simultaneously at 64 
positions, at various Doppler frequencies. The required transmitter activating signal is only 
30mVrms (on certain piezocrystals) , where the response of sound picked after propagation 
and detection is around SOOmVrms on the acquisition card. 

The time sequence of the measurement process is implemented with two National Instru- 
ments pci6071 cards installed on a PC running on Labview, and it goes as follows: A software 
key initiates once the sequence, in which a rectangular pulse is sent repeatedly and activates 
a switch that connects the main signal to the acoustic emitter. The pulse repetition rate is 
505//S and its duration IS/is. The main signal is sinusoidal with implied frequency between 
220KHz and 6MHz. The emitted sinusoidal signal must remain coherent between pulses 
and this will not be provided in burst mode of some function generators. The rising edge of 
the pulse triggers additional three clocks (two event clocks are available in each card). 1) A 
sample pulse is sent to engage the sample and hold components in the 64 lock-in amplifiers. 
The pulse width is 8/xs, which defines the integration time; delay time of 204/xs is chosen 
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according to the acoustic propagation time with additional 5//s to allow for geometrical de- 
lays, time shifts and sound from secondary sources to arrive on the detectors. The end of 
the sample pulse triggers the acquisition scan by the cards. 2) A gain variation pulse is sent 
to the preamplifier during the sample pulse and l/is before in order to increase the gain to 
a preset maximum value. Such step is taken to decrease the chance for a preamplifier to 
saturate and hang during the sample time in response to sudden strong noise (e.g. reflection 
of side pointing beam from the wall) . 3) A clear pulse is sent to reset the sample and hold 
components by opening a path to discharge their capacitors. The discharge starts after that 
sufficient time is allowed for the acquisition cards to scan the channels (scan rate of SOOKHz 
was used, and a waiting time was set assuming the cards work in sequence). The discharge 
time was set to 40/xs (twenty times the characteristic discharge period). 
In order to avoid errors by some arbitrary offset in the sampled signal a strobe mechanism 
was utilized. The frequency of the two reference signals, which enter the product detectors 
in the lock-in amplifiers, was shifted by lOOHz in relation to the main signal (coherence 
remains since the three function generators were connected on the same lOMHz time-base). 
The sequence of pulses is not synchronized with any of the function generators hence an 
initial arbitrary phase has to be resolved. In addition, since Labview on a PC is not a 
real-time software it is not promised that the two cards begin acquisition in the same pulse. 
The arbitrary phases are resolved by sacrificing two channels in each card. Thus, 62 lock-in 
amplifiers are used for acoustic detectors and another two are connected to the main signal 
(attenuated). During analysis the phase of the latter is subtracted from the former depend- 
ing on which card the lock-in amplifier is sampled. The direct hashing of the acquired data 
that provides the amplitude and phase of sound is based on summing the signal around the 
peak in Fourier domain. The center frequency is lOOHz as expected and the characteristic 
width is 50Hz. A summing window of about lOOHz is used, which is enough to include 
Doppler shifts by the time changing flow. Note that a zero Doppler shift is expected for 
steady state flow (although locally the wavelength is changed by the velocity of the flow). 
The given pulse repetition rate determines that only frequency shifts larger than 1800Hz 
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could alias with the main frequency signal. 

Significant part of the efforts in building the probe were devoted to design a suitable acoustic 
emitter. The important parameters are the length of the transducer in the axis parallel to 
the measurement plane, and the smoothness of the phase and amplitude. A long transducer 
is necessary to produce a beam of parallel rays, to minimize the effect of side lobes, and to 
reduce the edge diffraction on the detector. The transducer lengths are varied between 3mm 
and 145mm. The transducer width is around 10mm, close to the width of the detector and 
determines the thickness of a layer, in which the fiow is effectively averaged. 
A composite piezoelectric material was used in order to decrease lamb waves distortion. How- 
ever, by using piezo-composite in the largest sizes a smooth phase still could be achieved 
but not the amplitude. The following types of transducers were built: 
(Tl) This type was used for the transducers of 3.5mm and 10mm lengths. A piezoelectric 
plate (SPC-Y02) was covered by 4mm flat lens of Perspex and matched by conductive epoxy 
(Circuitworks) that was used also for wiring. The backing was air. 

(T2) A piezoelectric plate (PZ34 from Ferroperm) was covered by a flat lens made of perspex 
with groves for the soldered wires. The aperture size was 35mm, and the backing was air. 
(T3) Two or three pieces of PZ34 plates were adjusted together side-by-side in a frame to 
make a length of 96mm and 145mm, respectively. The backing was air, and the cover was 
a thin layer of epoxy. Wires were soldered on the edges. 

(T4) A better transducer compared with T3 was built from polarized PVDF, of 125X10mm 

an active area, mounted on Balsa wood and covered by a thin, conductive layer of epoxy 
and a thin layer of Acrylic. 

In all types of the emitters the enclosure was made of black Delrin, with 30dB/cm absorp- 
tion. Types Tl and T4 were levelled with a round wall. The latter was made using a large 
opening made in the cell with house flUed of water and covered by 10/xm Mylar fllm. 
Calibration of the array detector response in each acoustic element to a plane sound wave 
was accomplished by moving a small emitter (19x6 mm^ at 2.5 Mhz from Automation In- 
dustries) slowly and precisely, paralleled to the array and taking measurements in steps of 
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1mm. Summation of the complex wave function over all steps yielded the response to an 
ideal large emitter. 

2. Particle Image Velocimetry (PIV) 

The standard part of the PIV system was built by Oxford Lasers. A double exposure 
1Mb CCD camera and a 15mJ pulse infrared laser providing light sheet are used to measure 
the velocity field of the flow in two dimensions, with particles as markers. The seeding 
particles are 20/xm Polyamide from Dantec, dissolved in a 50% mixture with a surfactant 
(PolyOxyethylene-Sorbitan monolaurate from Sigma). The plane of the velocity measure- 
ment is chosen to overlap the acoustic measurement plane. The analysis is based on cross 
correlation between images by Fast Fourier transformation (FFT) of divided regions in the 
frame (usually by 32X32 pixels). Analysis of each region provides one velocity vector, but 
the algorithm can take overlapping regions up to 75%, so that for regions of 32X32 pixels 
from a 1Mb pixels camera a mesh of 128X128 vectors is obtained. In the implementation, 
the pixel equivalent size was 80X210/xm, the laser pulse duration was set to 60/xs and the 
pulse separation was set to 2130//S. Thus, correlation between consecutive frames could 
fetch rehably velocities between 0.1-l.lm/s (in a worse case, using 32X32 pixels). There is 
a tradeoff between mesh resolution and velocity accuracy and limits. Our analysis expands 
the velocity limit by combining results of analysis with different area sizes of a correlation 
region. We limit our search to low velocity gradients and thus we can interpolate higher 
velocities in a fine mesh that reveals the structure of the lower velocities. 
The system was modified to use it in the narrow field of view available between the rotating 
plates. Since the design of the experimental set-up (in two disks configuration) does not 
allow the camera to look right above the light sheet plane, we were accounted with three 
problems, (i) The view angle should minimize pickup reflections from the wall, (ii) The 
converging effect of the circular wall of the water filled cell as well as reduced optical clarity 
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can be corrected by using windows made of flat pockets of Perspex flUed with water outside 
the ceU. For reducing reflections separate windows for the laser sheet and the camera were 
found preferable, (iii) Focusing on a plane in a perspective view should use the Scheimpflug 
condition, in a way known as a perspective control. The perspective control requires break- 
ing the alignment between the camera and the lens so that the plane of the CCD (image 
plane) meets the principle plane of the lens at the same level, where the object plane meets 
with the other principle plane of the lens. The required tilting of the lens comes on the 
expense of reduced light intensity entering the CCD camera. In order to increase the light 
intensity we used two lenses. The tilted lens chosen was Mamiya RZ67 75mm/4.5 short 
barrel lens, capable for the perspective control. On the digital camera, we flxed a micro 
Nikkor-Nikon 50mm/l:1.4 lens through a F-C mount, which was shorten down to 11mm. 
This configuration allows us to focus on almost 15cm of the plane from 30cm distance, sub- 
ject to setting obtuse viewing angle. 



Following a derivation of an equation for sound scattering and refraction by a flow mostly 
due to Kraichnan [3] and Lighthill [21] one gets in the flrst approximation of a small Mach 
number, M — v/c <^ 1, the following equation for the sound scattering caused by the 
velocity fleld alone: 



Here ^ is the complex wave function that represents the sound pressure oscillations 
generating at frequency u by the emitter, ko = tu/c is the sound wave number, v and v' are 
the flow velocity and the velocity oscillations due to sound propagation, respectively, and po 
is the uniform fluid density. Eq.(l) was obtained in the approximation where low frequency 



III. THEORETICAL BACKGROUND OF THE DATA ANALYSIS 



A. Sound scattering by a finite width beam. 




(1) 
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sound generated by the flow and the sound attenuation due to viscosity were disregarded. 
Further, we assume 2D symmetry for simphcity, and due to the fact that the sound frequency 
is much larger than any frequency in the flow a quasi-stationary case is considered. The 
influence of, e.g., moving vortices in the flow shows up in a frequency shift (Doppler shift) 
relatively to the incident wave frequency. The former usually is smaller or comparable with 
characteristic frequency of fluid flow and is not considered here. In order to evaluate Eq.(l), 

— * 

we choose the axis xi along a local incident wavefront direction, thus || xi , and expand 
the source term. There is no influence of the curvature in the wavefront direction on the 
derivatives within the first order in the Mach number approximation [22]. Then taking into 
account that the wave length is the smallest scale in the flow one gets for an incompressible 
flow the following equation (in a plane wave approximation) 

A^ + kl^ ^[2klv-x-2iko^^^^^]^/c, (2) 

where Xi is replaced by x and is used to define the local ray direction. 
One can define the sound scattering field as "^scat — — "^rest, where "if rest is the complex 
wave function in the absence of a flow that satisfies the following wave equation without a 
source term: 

Arrest + kl^rest = 0. (3) 

By applying further the Green function method the solution of the scattering problem 
is obtained for a two-dimensional geometry in the following form (in ko\r — r'\ ^ 1 approx- 
imation) : 

dvx sGxp {iko\f — r'\ 



*scat = - exp (i7r/4) / (fr'{klv:c - iko^)^^=^L 
c J ox ^2T:ko\ 



-'^rest- (4) 



|r — r I 

In a far-field the solution can be drastically simplified and reduced to the analytical 
relation between the scattering field and the Fourier transform either of the velocity or the 
vorticity fields. Indeed, in the far-field \r — r'\ ^ \r\ — f ■ r', and the integral in Eq.(4) 
becomes the Fourier transform of the velocity and the velocity gradient fields. The variable 
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in Fourier domain is the scattering wave vector ks — ko(f — x) and = 2A;osin^/2, where 
f is the unit vector from the center of the scattering region toward the detector. The Fourier 
transform in 2D is defined as 



Taking into account that Fk^{^} — i{ks • x)Fk^{vx} — iko{cos9 — l)Fk^{vx}, one obtains 
in the far-field hmit 

1 (27rA;o)^exp(i7r/4) . ~ . , , 

*,eat = ^ COS 9 Fk.iVx"^ rest}, 5 

where cos 6' = f ■ x and '^rest = ^rest exp (— i/co;'). We present in eq.(5) the scattering 
field as a convolution of the Fourier transforms of the velocity field and the known beam 
function, Ffc,{^rest}- Only in an ideal case of an infinite planar wave front, when ^rest — 1; 
there exists a simple relation between the Fourier transforms of the velocity and the vorticity 
fields: 

= %^i^..{(V X v%} = ^ cot m)FkAi^ X vU, (6) 

which leads to a direct linear relation between the sound scattering field and the Fourier 
transform of the vorticity field as first pointed out by Lund and Rojas [5]. 

This general approach can be applied in a case of a single rigid body rotation vortex 
with a core azimuthal velocity distribution as 

vt{r) = fir, for r < ro, (7) 

and outside the core up to the cell wall the azimuthal velocity component decays as 

^«(^) = ^ — ^(-^ + ^c/^)' for ro<r < r^, (8) 
where ft is the angular speed, Tc and Tq are the cell and the core radii, respectively. 
In the axisymmetric case one gets the following result for the Fourier transform of the x 
velocity component, Vx'- 

FkA^x] = 77- COS {9/2)Q-——^^[Ji{ksro)rc - Ji{ksrc)ro\, (9) 
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where Ji{x) is the Bessel function, which was introduced via Fk^{vx} — 
~^2tk^^ / r'vt{r')Ji{ksr')dr'. The expression (9) exhibits even for the unbounded beam width 
and detector length two maxima in the scattering signal at some value of kg due to a final 
size of the flow region, Tc- It is clear from Eq.(9) that in a general case the location of the 
maxima depends on two parameters, ro and Tc- So in the case when tq/tc <S 1, the peaks 
are located at kgTc — 1.1 Itt. This value alters by change of Tq. So the numerical calculations 
based on Eq.(9) give the following results (for Tc — 145 mm): at ro = 6mm the peak position 
is at ksTc — I.IItt, at ro = 20mm the value is kgVc — 1.097r, at ro — 50mm the expected value 
is ksTc — 1.027r, but at ro = 100mm the location value is already shifted to kgVc — O.STtt. 
The peak height, H, also provides useful information about the flow. Indeed, in the case 
of ro/rc <S 1 and of the inflnite beam width, H is proportional to Q, to a core circulation, 
r = and independent of fco; or sound frequency. The dependence on T and indepen- 

dence of /co are rather non-trivial results taking into account rather complicated functional 
dependence of Fk^{vx} on the parameters ro, rc, and k^. 

In the limit rc — > oo the peak locations approaches zero angle asymptotically but the value 
of the amplitude of the Fourier transform of the velocity at ^ = remains zero. In the limit 
of kgTc <^ 1 and at ro <S Tc one can derive an asymptotic expression for the Bessel functions 
in Eq.(9) and get Fk^{vx)}\ksrc^o = —^r^rlko sin 9/2. From this expression it is easy to see 
that the amplitude at small angles increases proportionally to fl, F, and kg, and approaches 
zero at ^ — > 0. 

For a finite width beam (or emitter) or a finite width detector the Fourier transform of 
the scattering field wave function from Eq.(5) can be expressed as a convolution ^scat(^) ~ 
Fksi'^x} * Fk^i^rest}- In ^ case of the finite width, d, (either of a beam or a detector what- 
ever is the smallest) let us consider first the Gaussian beam model. The beam function is 
defined as ^rest — exp (— y^/[2(ci/2)^]). For d/vc ^ 1 the behaviour of the peak height here 
is similar to the infinite width beam considered above. 

The location of the peaks in a finite width Gaussian beam depends on the beam width: 
a,t d — 50mm the peak position is at kgd ~ 0.947r (at Vo/d <^ 1). When ro increases this 
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relation is changed: at ro = 10mm the peaks are located at kgd ~ 0.927r, and at ro = 50mm 
they are found at kgd ~ O.TItt. The dependence of the peak heights on the parameters for 
a finite width beam is more elaborate, since here three characteristic lengths exist in the 
problem: ro, Tc, and d. In spite of this fact the approximate scaling law was found from 
numerical calculations. It follows that H is independent of fco ^-nd proportional to F with 
the scaling function g[x) (where x = ro/d) that is different for two regions of small and large 
beam width: (a) g{x) — 125x^ '^ for x > 1, and (b) g{x) = 66 + 66x^ for x < 1. 
In the case of a finite rectangular beam (or a finite detector whatever is the smallest) the 
peak location is defined by kgd as in a finite Gaussian beam considered above. At ro/c? <C 1 
one gets kgd ~ 1.527r (e.g. for d ~ 50mm). For larger values of Tq the peak location is shifted 
towards smaller values. So, at ro — 10mm the estimated peak position is at kgd ~ 1.487r, 
and at ro = 20 and 50 mm one gets I.SQtt and 1.347r, respectively. The peak height similar 
to the case of the Gaussian finite width beam, is independent of ko and proportional to F 
with the scaling function g{x) that is slightly different from the former case. It is determined 
for the range of arguments values 0.01 < x < 4 ai d > 30mm and Tc = 145mm and has the 
following form: (a) g{x) = 205x^ '^ for x > 0.5, and (b) g{x) = 63 + lOSx'^ for x < 0.5. 
In order to demonstrate the influence of a finite beam width on the scattering signal from a 
single vortex with a core radius ro, we present in Fig. 3 the results of numerical simulations 
for the structure functions of the velocity field for the vortex of ro = 5mm and three beam 
widths of 20, 60, and 120 mm (but all of them smaller than the cell size, Vc — 145mm) 
compared with an infinite beam width for the same cell size. Only for d = 120 mm the 
wave number of the structure function peak location becomes close to the peak location for 
the infinite width beam (see Fig. 3a). The wave number of the peak location in the infinite 
width beam case is defined by the cell size (see explanation above after Eq.(9)). One can 
see that in spite of the fact that o? ^ ro, the finite width of the beam drastically alters the 
resulting scattering signal. Thus in order to get correct results on the structure functions of 
the velocity and vorticity fields one should use either a beam width exceeding the flow size, 
rc, or to perform a deconvolution based on a provided constrain, which fills out the missing 
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information about the velocity field outside the beam extent. 



B. Alternative approach to find velocity and circulation: phase shift slope in a 

forward scattering. 

While the study of scattering data requires to extract from sound signals both the am- 
plitude and the phase variations, one can find information about the velocity field and the 
circulation from the phase shift. We consider two cases in a phase shift analysis. In the first 
case of a large emitter, rays remain parallel in propagation, and the phase shift as function 
of y (assuming it can be measured avoiding interferences between the paths) can be defined 
in the hmit M < 1 as [22]: 

0(y) = (A;o/c) J v,{x' ,y)dx' . (10) 

Lindsay has found [2] that the phase shift induced by a point vortex has a linear depen- 
dence on an angle. This result can be obtained if we consider rays passing outside the core, 
and assume a large emitter at x — > — oo and velocity profile 1/r, namely: 



rji-Z _|_ yZ 



Phase shift is measured on a screen aA, x — re, with 6 = tan~^(^), then the slope can be 
found, outside the singularity at ^ = 0, with the following constant value: 



r ^dx = -^Qrl (12) 

J-oo C C 



89 de J-oo c 

Our idea is to look in particularly on rays passing near the center of the vortex at 6* = 
and try to get information on the fiow from the phase shift slope. For an infinite source at 
X — —Tc and detectors &i x — Tc- 

•90, /core dv^ 



89 c J-r-c dy 

Thus, substituting the velocity field projected on x via Eqs.(7,8), one obtains the phase shift 
slope as 

— 0^0 = 41]ro 1 ■ . (14) 

o9 c rc + To 
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-dx at y ^ 0. (13) 



Comparing Eq. (14) and Eq.(12) we find that a rigid body rotation core induces a larger 
slope by an order Tc/ro compared with the exterior of the vortex. Thus for a point vortex 
at 6* = there is a step in the phase shift, having an infinite slope, revealed therefore as the 
Berry's phase [23]. 

A total change in the phase shift between two points on the acoustic screen can be 
calculated by a loop integral over the velocity field. Since a velocity along the detector and 
the emitter is zero, there is a simple relation between the circulation and the step in the 
phase shift followed from Eq. (10) [24]: 

Y^jvdl^ Y^{^<t>{2) - A0(i)), (15) 

where A0(i) and A0(2) are the phase shift differences on two sides of the detector, respec- 
tively. 

In the second case we consider rays diverging out of a single point on a small emitter (and 
there is an axisymmetric single vortex). We can show that the phase shift slope, calculated 
from the phase shift due to a vortex at the cell center between a detector placed at x = Tc, 
and sound emitter located at (x, y) = (— rg, 0), can be written as 

koTc r= vt{r) 



o9 c 

The angle 9 is measured between the radius from the center of the vortex and the beam 
direction. The phase difference used to derive Eq. (16) can be written as [2] 



r^dr. (16) 
Jo r 



A(i){e) = ko r {n - l)ds', 

J—rc 



where n — ^1 - 2sv/c is the refraction index in the first approximation in the Mach number, 
s is the unit vector in the ray direction, and s' is the coordinate in the ray direction. Then 
from simple geometrical considerations one gets Eq. (16). Substituting the expression for 
the velocity profile Eqs.(7,8) into Eq. (16) we obtain: 

— 9^0 = 2flro{l ■ (17) 

09 c rc + To 

that is equal exactly a half of the result for a large emitter. 
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C. Construction of a far field wave function from the near field measurements. 



Probably, the most severe limitation of the theory to be applicable to an experiment is 
the far-field approximation. In this respect the relevant question is how to extrapolate sound 
scattering measurements obtained at a distance not far away from a scattering region into 
a far-field region to be compared with the theory and to reliably extract information about 
the vorticity structure function. Our validation of the far-field scattering result is based on 
comparison with calculations, which use the velocity field measured by PIV technique. 

A far-field construction of the scattering field from the acoustic field (either pressure or 
scattering wave function) given at a certain plane in 3D case and at a certain line in 2D case, 
is based on the mathematical description of a Huygens principle typically used to describe 
radiation from a curved surface. This problem is similar to the problem of diffraction by a 
thin screen of finite dimensions (or for a finite dimension wave beam) [25]. Thus our goal 
is to consider the propagation of a wave into unbounded source- free half-space x > 0, when 
certain conditions in the initial plane x — are specified. Then the scattering field in 2D at 
the location x >0 can be defined from the following Rayleigh-Sommerfeld integral [25] : 

■qf{x,y) = - 1^^^ ^{0,y') — [G{x,y;0,y')]dy', (18) 

where G{x, y; x', y') = G{\f—r'\) = ^Hl^'^\k\f—r'\) is the Green function of the Helmholtz 
equation in two dimensions(2D), r = {x,y), and Hq^^ is the Hankel function of zero order. 
Equation(18) follows from the integral Helmholtz theorem, where use of the Green's function 
is made. S+ and E_ denote the two sides of the diffraction screen {x > 0) and {x < 0), 
respectively. The presence or absence of the screen is not essential to the derivation of eq.(18) 
but in the absence of the screen one does not need to apply the Kirchhoff approximation. To 
get a more simplified expression in a far-field region the asymptotic expression for the Hankel 
function can be used Hq^\z)z^i ^ y^expi(^ — 7r/4). Then the far-field construction of 
the scattered sound wave function can be calculated in the limit of A;o|r — r'| ^ 1 and 
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ko(y - y'Y/r < 1 as 

y)ff-scat = koVi I dy' ' exp (^M^^^)^(o, y')d-.cat, (19) 

•> ^27rfco(r/ - Td) ^l^/ ~ ^d) 

where ra and r/ are the distances measured from the cell center till the detector and the 
far- field region, respectively, ^(0, y%_scat ^-nd '^{x,y)ff-scat ^-re the scattering wave func- 
tions at the detector (as measured) and at the far-field, respectively. Since at the array 
detector there are just 64 elements, a sum instead of the integral is used to calculate the 
wave function in the far-field. 

The method presented is found advantageous over other inverse solution methods. In mea- 
surements as well as in numerical simulations it is also desirable to apply tapering of the 
beam edge at the emitter to avoid an oscillating edge diffraction pattern. 



IV. EXPERIMENTAL RESULTS. 
A. Method of data analysis 

In order to get information about the scattering sound amplitude one needs to find a 
phase difference between an incident and a scattered signals. Our way to find the scattering 
wave is simply to measure and subtract the complex sound signals with a fiow and at a rest. 
At high enough frequencies we can concern ourself only with the main change in the wave, 
which is the phase shift due to change in time of flight by the flow. The phase shift is found 
from the phase information as 4> — — arg(^^^). 

The amplitude of scattering wave on the array detector is determined as ^4 = {"^fiow — 
"ifrestl- Both wave functions, "^fiow and "if rest are the central frequency component of a 
sound extracted from a bank of pulses entered to a Fourier type filter, as described above. 
Information obtained during a period when the bank of pulses is collected, is called a frame. 
In a steady flow we average the complex results between available frames. For the dynamic 
study of a fiow we compare the wave function in a fiow with the average on many frames 
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to avoid suspecting additional sources of fluctuations. Another step before substitution of 
"if rest into the last equations is to normalize it by a factor expi7. The phase correction 
7 is the average in time (between frames) of the values 7' that minimize the expression 
< l"^ flow — ^rest expi7'| >y (average on y refers to channels of the array). Such step is used 
to specifically regularize results of scattering from a single vortex and it also compensates for 
occasionally uncontrolled conditions that change the overall refraction (for example raising 
of seeding particles by the flow, when they are used). At the end, we obtain the scattering 
wave function on the transducer array plane as: 

scat flow -'^ rest- (20) 

The presented results of scattered and incident flelds are normalized so that the incident held 
on the emitter plane is of a unit amplitude. To calculate the incident fleld at a distance I 
from its plane we use the Huygens construction (see Eq.(18)) that is a propagation transform 
of the window function (with 5% apodization at the edges to reduce signal corruption due 
to noise) ^(0, y): 

'^rest{l,y) = ikoVi J dy'exp{iko^P + (y - y'^) ^^J-^ ^rest{0,y')- (21) 
Then ^scat can be calculated from Eq.(4). 



B. Characterization of the flow by PIV. 

We focus on a single vortex, created between co-rotating plates, and find its average 
properties in time. Since the vortex is axisymmetric we can extract from PIV analysis 
the profile of an azimuthal velocity versus the radius. The center of the vortex is found 
according to a minimum in the orthogonal projection of a velocity vector on x and y axes 
(except for the smallest vortices produced by 25mm rods, where we designated the center 
by eye recognition). The profile is built by averaging on 64 PIV maps, such as shown in 
Fig.4a. The profiles for the three inspected fiows are shown in Fig.4b (co-rotating 25mm 
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rods), Fig.5 (co- rotating 40mm plates), and Fig. 6 (co-rotating 100mm plates). In general, 
we see that rotating plates with blades produce a core of a rigid body rotation, i.e., the core 
velocity increases linearly with the radius. Outside the core, the velocity decreases like in a 
flow between two cylinders. This solution (Eqs. (7,8)) was used to fit the measured velocity 
profile, and the flow parameters were extracted from the parameters of the fit. Thus we find 
from the fits ro = 6mm, ro = 14mm and tq = 42mm, respectively. In all cases generally 
the core rotation frequency is the same as of the motor. The PIV maps extend to about 
third of the cell due to geometrical limitations in a field of vision and a limited power of 
the laser (nominally 15mJ/pulse). Therefore, the velocity profile is extrapolated using the 
model of a flow between two cyhnders. The distance from the cell center to the detector's 
plane and the wall is Tc — 145mm. We used an average on 64 velocity maps of the flow for 
40mm plates (extrapolated between radii of 50mm and 145mm) to evaluate the integral of 
the orthogonal velocity projection on x-axis. The result can be compared to the ultrasound 
phase shift measured at 5. 5 MHz, as shown in Fig. 7. 

C. Velocity and circulation measurements via phase shift. 

As follows from Eq.(15) the measurement of the phase shift provides direct information 
about the circulation. This result can be obtained from the data in Fig. 7, where the spatial 
dependence of the phase shift obtained from the sound and PIV (here the velocity field 
measurements were converted into the phase shift via Eq. (10)) measurements are presented. 
We would like to emphasize here that all experiments presented and analyzed were realized 
in the limit M <C 1, so that the expression for the phase shift delivered in Eq.(lO) is 
applicable. It follows that for the maximum phase difference measured in Fig. 7 of about 
A(p — 0.7rad, the maximum circulation is about P — 45, OOOmm^/sec. As one can find 
from Fig. 7 the value of the maximum phase difference and correspondingly the maximum 
circulation coincide with those found from the PIV measurements within 1-2%. 
The phase shift patterns due to the flow based on averaging of 10 frames of 512 pulses for 
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various emitter sizes were examined. It appears remarkably for two types of flow considered 
(both a single vortex flow with a rigid core rotation, produced by 100mm plates and 25mm 
rods) that the slope of the phase shift around the center of the beam versus scattering angle 
(or a distance along the detector) is constant. We describe in detail the flow produced 
between two co-rotating 100mm plates, which provides stable phase shift plots. In Fig. 8 we 
demonstrate how a value of the phase shift slope d(j)/d6 are obtained: the limit is chosen 
to give a constant slope with the smallest linear regression error. According to Eq.(17) the 
phase shift slope provides direct information about the vorticity of the vortex, Co — 2Q. The 
results for emitters of different sizes and at various frequencies are shown in Figs.9a-d. The 
phase shift slopes for all frequencies and all emitters appear to be linearly dependent on 
the rotation speed, il. The data presented can be either scaled by Q or presented via the 
derivative on il, m = Then the data on the proportionality coefficient m as a 

function of frequency, calculated from the plots presented in Figs.9a-d, are summarized in 
Fig. 10. The proportionality coefficient, m, depends linearly on the frequency in the range 
between 0.5 and 5.5MHz, and the results are sharply separated in two groups: (i) small 
emitters (3.5, 10mm), and (ii) large emitters (35, 145mm) compared with the size of the 
detector array of L = 62mm. The rotation speed of the plates is used to define the flow 
parameter and is converted to the angular speed of the vortex core obtained from the fits of 
the PIV data. According to Eq.(17) in the theoretical section, we expect a point-hke emitter 
to have the frequency dependence of the proportionality coefficient, K — dm/df, on the 
plots of Fig. 10 as follows: 

^ - '^^W^ = = 0.0281(.ec)(MF.)-'. 

oiloj Vc + ro 

The experimental value of K for 3.5mm and 10mm emitters is K = 0.026 ± 0.001 
sec(MHz)~^. The value of K from the plot in Fig. 10 for 35mm and 145mm emitters is 
0.045 ±0.002 sec(MHz)~^. As explained in section IIIB the theoretically predicted value of K 
for a small emitter should be twice smaller than for a large one, i.e. the theoretically expected 
value for the large emitters is K — 0.056sec(MHz)~^ according to Eq. (14). However, 
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specifically for these measurements the large emitters blocked 12% of the total diameter 
of the exterior of the flow. Assuming effective reduction in the cell radius Tc the expected 
theoretical value should be corrected down to = 0.048 sec(MHz)~^, and it becomes rather 
well comparable with the experimental value presented above. 

D. Sound scattering data. 

Using Eq.(4) derived in the theoretical section we were able to calculate the scattered 
signal amplitude from a given velocity field assuming that the emitter wavefront is known 
exactly. We used a velocity profile of a flow produced by 100mm rotating plates extracted 
from PIV (Fig. 6). The emitter of 35mm length was used and modelled for a flat amplitude 
(a plane wave) as the input, to obtain the result in Fig.lla. The measured incident and 
scattered wave amplitudes at the frequency 4.0MHz are compared to the calculated ones 
(they are scaled to unity value of the wave amplitude at the emitter exit) (see Fig.lla,b). In 
the calculation of the sound scattering field, a mesh of a half wavelength resolution was used 
in a zone limited to twice of the emitter size and the velocity field was interpolated based 
on the velocity profiles values. The wavefront was calculated at various cross-sections in the 
cell using Eq.(21). The measurement was performed at 5.5Mhz collecting 20 x 512 pulses 
in a flow produced by the 100 mm disk at the rotation speed Q — 22 rad/sec. In Fig.l2a,b 
similar data of the incident and scattering sound fields from 125mm emitter at the frequency 
5.5MHz are shown (100 mm disks, fl — 35 rad/sec). As in the previous case an agreement 
with the theoretical calculations is rather good. An evidence of a diffraction pattern (side- 
lobe) from the emitter of 3.5mm length at the frequency of 5.5MHz in still water (regarded 
as the incident wave), and in the scattering signal due to vortex flow between the rods of 
25mm diameter at Q = 71rad/sec, is clearly seen in Fig.l3a,b. The direction of the initial 
sound beam was slightly tilted relatively to the direction to the cell center, so that the center 
of the beam is shifted from the center of the flow. This fact was not regarded in the fitted 
curve that results in a discrepancy of the measured and calculated sound amplitudes. This 
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result clearly demonstrates that the scattering signal can be easily buried in the incident 
signal coming from the side-lobes at sufficiently large angles of detection. 

E. Far field construction of sound scattering field and its relation to velocity and 

vorticity of the scattering region. 

We use the Huygens projection of the near-field scattering signal detected at the re- 
ceiver array plane, into a synthetic far-field plane using the Rayleigh-Sommerfeld integral 
via Eq.(19). We choose the far-field plane at a distance of r// = 250m, much larger than 
(Pko/STT = 3m. The signal is tapered 2mm on each edge before the projection to suppress 
numerical instability of the diffraction pattern. By this procedure one obtains full informa- 
tion of an amplitude and a phase of the scattered signal in the far-field. However, practically 
due to finite spatial resolution the phase information becomes corrupted. So the next step 
in the far-field construction is to replace the phase variation as a function of an angle (or 
a wave number) by a function toggled at each minimum point by a ±7r/2. In such way we 
reconstructed the phase field using an observation made in the simulations that at every 
minimum point of the correct amplitude curve the sign of the field should be inverted. In 
order to avoid some spurious minimum points the number of digitizing points in integration 
was increased, and the result was compared against finer digitization mesh. 
According to Eq.(5) the sound scattering field in the far- field is proportional to the Fourier 
transform Fk^{vx^rest}, in such a way plots of the modulus of the two-dimensional 
Fourier transform, \Ff;^{va;^rest}\: in the entire cell were obtained in Figs. 14, 15. The data 
projected into the far-field are compared with calculations of the scattering field based on 
the velocity field measurements by PIV. The results are presented in a view angle observed 
from the center of the cell through the receiver as an aperture. This angle, 9, is related 
to the scattering wave vector, kg, via formulas: kgx = kQ{cos9 — 1), kgy = kosinO, and 
ks = 2kosm9/2. 

We studied properties of the modulus of the Fourier transform, \Fk^{vx^rest}\: different 
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sound frequencies and different rotation speeds. It was revealed tfiat tfie angular location 
of the main peaks in \Fk^{vx'^rest}\ is proportional to the sound wavelength (Fig. 16). The 
experimental value of the slope of the plot is 2 sin (6*^/2) (/cqc?) = ksd = (1.51 ± 0.11) • tt, 
and is found to be in a good agreement with the value of the slope obtained from numerical 
simulations kgd — 1.52 • tt (see discussion in Sec.IIIA). It was also found that the scaled 
value of the peak height P = Hd~^'^g{rQ/ d){14:b/rcf''^ of \Fk^{vx^rest}\ is proportional to 
the circulation F (Fig. 17). Thus, four different sets of the data for three different vortices 
and two beam sizes can be scaled down using the functional dependencies on the beam 
width, (f-^, and the reduced vortex core size, g{rQ/d), taken from our numerical calculations 
(see discussion in Sec.IIIA). As seen from the plot scahng does not work great for one set of 
the data, possibly, due to flow imperfection. On the other hand, the proportionality of the 
scaled peak height to F is not so obvious as one can decide from the first sight, since the finite 
beam width described by the beam function ^rest, alters significantly the function presented 
in Eq.(9). We also found from the measurements that the peak height, H, is independent 
of the sound frequency and depends strongly (about when d < 2ro) on the beam width. 
When d » 2ro the peak height value, H, is proportional to Tq in a full agreement with the 
simulations. It is obvious that the finite size of either a beam (an emitter) or a detector 
(whatever the smallest) limits our knowledge about the scattering field and, therefore, also 
about the Fourier transform of the velocity field at the wave numbers (or angles) smaller 
than that corresponding to the peak location.. 

The projected scattering to a far-field, Fk^{vx'^rest}, does not have a direct relation to the 
structure function of a vorticity, Ffe^{(V x v)z}- In order to get the latter, one should extract 
information about a 2D Fourier transform of the velocity field Vx from the measurements 
of Fk^{vx'^rest}- Howcver, this transformation is a singular one, so a choice of constrains 
is made roughly related to the continuity and irrotaional character of the fiow outside of 
the acoustic beam extent, where ^rest decays to zero. Particularly for the plane beam, 
a ID Fourier transformation is relevant, since in the propagation direction, x, Fk^i"^ rest} 
can be presented by a delta-function in the wave number domain. We performed a back- 
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ward and forward ID Fourier transformation of Fk^{vx'^rest} values on a set of points in 
9 coordinate (which span about 90deg), such that the values are not modified initially. 
Next we added the filtering effect that removes "if rest by extrapolating the intermediate re- 
sult Vx{y) = 27r / dhsyBxp {iksyy)F {v^"^ rest} ^ if it were one-dimensional projection of the 
velocity, vl.^{y), using our flow model for extrapolation (see Fig. 18). The phase in the ex- 
trapolated region of Vx{y) should be adjusted having continuity at y = ±d/2 and a tt phase 
difference at the tails. Such construction provides the required velocity field and correspond- 
ingly its Fourier transform in the wave number domain, F{vx} = J dy exp {—iksyy)vx{y). 
Now using Eq.(6) one can relate in such a way the obtained Fk^{vx} to the required structure 
function of the vorticity -Ffe^{(V x v)^} (see Fig. 19a). With a help of the Wiener-Khinchin 
theorem one can also obtain an azimuthally averaged point by point correlation function of 
the vorticity, which in the case of an isotropic flow or axial symmetry gives 



The corresponding correlation function of the vorticity is shown in Fig. 19b. Since the weight 
of errors in estimation of F{vx} is more pronounce at high scattering angles, the integral in 
Eq.(22) should be cut in a tight range. We used a Gaussian attenuation factor (with kg — 
0.7mm~^ characteristic value) on the far- field scattering values to filter out high frequency 
fluctuations. 

It is obvious that the vorticity structure function restored in such a way contained additional 
information reflecting our guess (or our PIV results). However, this additional information 
is mostly relevant to a low wave number range of the vorticity structure function, i.e. to the 
core of the Fourier transform of the vorticity. Therefore, at wave numbers larger than the 
core wave number value, the vorticity structure function provides information based on the 
original experimental data. 




(22) 
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F. Dynamics of a vortex measured by sound scattering. 

The sound scattering technique can be used to study vortex dynamics. It can be studied 
by a phase dynamics approach to get temporal variation of a vortex location, vortex radius, 
and vortex circulation similar to what was done in Ref. [16,17]. However, we used a different 
approach of sound scattering and compare the results on a vortex precession with those 
obtained by PIV. The vortex position is found from the minimum point of the scattering 
amplitude pattern. The periodicity in the vortex motion is characterized by the variation of 
the peak heights in the far-field scattering pattern. The setup is used with the upper plate 
of 280mm in diameter at the angular speed Q = 30rad/sec, and the lower plate is absent. 
The rate of PIV tracking is one map per 0.266 seconds. The rate of ultrasound tracking 
based on collection of 32 pulses per frame is one wave function plot per 0.016 seconds (sound 
frequency is 2.5MHz, emitter is 96mm long). In the case of the rotating upper plate, a pe- 
riodic precession of the vortex is found by both techniques: PIV and the sound scattering 
(see Fig. 20). To compare two sets of the data the time correlation functions of the vor- 
tex locations were produced from both sets. A good quantitative agreement between the 
dynamic results of PIV and ultrasound measurements was found (see Fig.21)with a period 
of O.QOsec. Besides, the peak height of the Fourier transform of the velocity projected into 
a far-field is found to be twice periodic compared with the vortex position (see Fig. 22a). 
It is clearly seen in the auto-correlation function presentation of the data in Fig. 22b. This 
effect occurs due to increase or decrease in the integral of Vx over the beam area (detector 
view), when the vortex core is shifted sideways, since for small kg \Fk^{vx'^rest}\ contains a 
contribution from the integral of v^, and double periodicity shows up due to absolute value 
of the function. The precession frequency as a function of the rotation frequency obtained 
from the auto-correlation functions similar to that shown in Fig. 22b, is presented in Fig. 23. 
We were not able to get such features from the PIV measurements due to restrictions in 
accuracy. 
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V. CONCLUSIONS. 



We demonstrated that our system provides rather rehable information about both the 
phase and the amphtude of the sound scattering signal by acquiring simultaneously 64 
channels of the detector array. The spatial and temporal information on the phase of the 
scattered signal allows us to get values of circulation, vorticity, vortex location, and vortex 
core radius. We verified quantitatively the theoretical value for the slope of the proportion- 
ality coefficient, K, for a single rigid body rotation vortex in a finite size cell. This method 
is rather comparable with the acoustic time-reversal mirror (TRM) method [15-17]. Instead 
of vorticity (or phase difference) amplification due to number of crossing of the flow in the 
TRM method, in our method the signal-to-noise ratio is amplified by averaging complex 
wave functions over many pulses. So the comparable resolution in the phase difference mea- 
surements can be achieved. 

At the same time in our approach we can also use the amplitude of the scattering signal to 
characterize the flow. The existing theories of the sound-flow interaction provide a relation 
between the scattering signal and the structure functions of velocity and vorticity flelds of 
a flow. However, strong limitations imposed by the theories make an application of the 
sound scattering technique for flow characterization rather restricted. The theories consider 
the scattering signal produced by a sound emitter and obtained by a receiver of an inflnite 
length in a far-field limit. 

In the paper presented we studied both theoretically and experimentally a possibility to 
obtain reliable information about velocity and vorticity fields of a single either stationary or 
time-dependent vortex fiow by the acoustic scattering technique with a finite width sound 
beam of the order of the vortex size and a finite size receiver taking a scattered signal rather 
close to the scattering region. The experimental results and analysis are compared with 
the PIV measurements taken simultaneously on the same vortex. From the theoretical side 
the main step in our analysis is the use of the Huygens construction to obtain the acoustic 
scattering signal in a far field from the known sound distribution in any intermediate plane. 
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This reconstruction allows us to overcome rather difficult experimental problem to obtain 
directly the scattering signal in a far-field. Such procedure also helps us to reafize that the 
only empirically relevant calculation of scattering is that of a confined fiow, since an emitter 
and a detector define the fiow perimeter. 

Another theoretical suggestion is the use of the beam function ^rest, which describes the 
finite width beam and allows us to use the same formalism as for an infinite width beam. 
The analysis of the experimental data based on the revised theory of sound scattering shows 
rather good agreement between the sound scattering and PIV measurements, i.e. the sound 
scattering signal obtained from the velocity field measured by PIV, coincides rather well with 
the scattering signal obtained in the ultrasound experiment. The same can be said about 
the Fourier transform of the vorticity obtained from the scattering via the construction sug- 
gested and that reproduced from the PIV measurements. We also show that the peak value 
of the far field scattering signal is proportional to the circulation or the angular speed of the 
core of a single vortex, and the angle distance between the peaks is inversely proportional 
to the sound frequency. Further natural step is to use this approach to a turbulent fiow and 
to study the vorticity structure function in a turbulent fiow. 

This work is partially supported by Israel Science Foundation grant, by Binational US- 
Israel Foundation grant, and by the Minerva Center for Nonlinear Physics of Complex 
Systems. 
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FIGURES 

FIG. 1. Experimental set-up: (1) rotating plates, (2) light sheet and scattering plane, (3) 
PIV camera with Schleimpflug optic arrangement, (4) acoustic emitter, (5) main signal pulse, (6) 
transmit switch, (7) detector array, (8) 64 preamplifiers, (9) 64 lock-in amplifiers, (10) reference 
signal, (11) sample and hold trigger, (12) acquisition cards (128 channels). 

FIG. 2. Geometry of cell, fiow and sound scattering. 

FIG. 3. Numerical simulations of the Fourier transform of the velocity field of a single rigid 
body rotation vortex flow vq = 5mm for a finite width beam, d, compared to an infinite width 
beam: (a) d = 120mm, (b) d = 60mm, (c) d = 20mm. Solid line-a finite width beam, dash-dotted 
line-infinite width beam. 

FIG. 4. (a) Velocity field snapshot obtained by PIV for a single vortex flow created by 
co-rotating rods of 25mm diameter; (b) Azimuthal velocity as a function of a cell radius extracted 
from PIV measurements for a flow between rotating rods of 25mm diameter at Q = 71rad/sec: 
black dots-experimental points, solid line-the model with fitting parameters: the angular speed 
Q = 71rad/sec and the core radius Vq = 6mm. 

FIG. 5. Similar to the plot in Fig.4b but for a fiow between rotating plates of 40mm diameter 
and n = 42rad/sec. Solid line presents the fitted model with the core radius ro = 14mm 

FIG. 6. The same for a flow between rotating plates of 100mm diameter at O = 23rad/sec. 
Solid line presents the fitted model with the core radius ro = 42mm. 

FIG. 7. Spatial dependence of the phase shift by a flow between rotating plates of 40mm 
diameter at = 42rad/sec: stars-from sound measurements at 5.5MHz, circles-from PIV. 

FIG. 8. Phase shift as a function of the scattering angle. The solid line is the fit, by which the 
the phase shift slope is defined. The data are for a fiow between the plates of 100mm diameter 
and J7 = 23rad/sec with an emitter of 35mm length at / = 4MHz. 
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FIG. 9. (a) Phase shift slope scaled by sound frequency, {d<p/d6){l/ f), as a function 
of angular velocity for a flow between the plates of 100mm diameter and with emitter 
of 35mm length; (b,c,d) Phase shift slope {d(p/d9) as a function of angular velocity for 
a flow between the plates of 100mm diameter and with emitters the of 143mm length, 
10mm length, and 3.5mm length, respectively, at various frequencies: squares-5.5MHz, di- 
amonds-4.5MHz, dots-4.0MHz, left-triangles-3.5MHz, up-triangles-2.5MHz, hexagons- 2.0MHz, 
stars-1.5MHz, down-triangles-l.OMHz, right-triangles-0.5MHz. 

FIG. 10. The derivative of the phase shift slope on the rotation speed, m, as a function of 
sound frequency for a flow between the plates of 100mm diameter for different size of emitters: 
squares-3.5mm, dots-lOmm, stars-35mm, triangles- 143mm. Solid lines are the linear fits. 

FIG. 11. Modulus of incident (a) and scattering(b) sound amplitude spatial distributions ob- 
tained by the receiver (black dots) are presented. Solid lines are the calculations. The flow is 
produced between two rotating plates of 100mm diameter, the emitter of 35mm length at 4MHz 
is used. 

FIG. 12. The same as on the plots in Fig.lla,b but for the emitter of 125mm length at 5.5MHz. 

FIG. 13. Evidence of side-lobs (diffraction pattern) obtained with emitter of 3.5mm length at 
5.5MHz: (a) incident, and (b) scattering sound amplitudes. 

FIG. 14. Modulus of the structure function of the velocity and the beam function projected 
into a far field for a flow between the plates of 40 mm diameter and with the emitter of 125mm 
length. Dots are the experimental data from sound scattering, solid line is calculation based on 
PIV measurements with the flow parameters = 42rad/sec and ro = 14mm obtained from the fit, 
and 60 mm characteristic beam size is assumed due to the aperture. 

FIG. 15. The same as on the plot in Fig. 14 but for the emitter of 35mm length and the flow 
parameters = 23rad/sec and ro = 42mm obtained from the fit. 
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FIG. 16. Angular location of the structure function peak, 2sin(^p/2), as a function of the 
reduced sound wave length, {kod)~^. 

FIG. 17. Value of the scaled structure function, \F{vx^}\giro/d)d~^-^{li5/rc)^-'^, peak, P, as 
a function of the circulation, F: stars- vortex flow with tq = 6mm created between two 25mm 
rods, 60 mm beam width, and emitter T4; dots- vortex flow with vq = 42mm created between two 
100mm plates, 60 mm beam width, and emitter T4; squares-vortex flow with tq = 14mm created 
between two 40 mm plates, 60 mm beam width, and emitter T4; left-triangles-vortex flow with 
ro = 42mm created between two 100 mm plates, 35 mm beam width, and emitter T2. 

FIG. 18. Reconstructed one-dimensional velocity field {vx)id according to the procedure de- 
scribed in the text for 60 mm beam width, for a flow between two plates of 40 mm diameter and 
Q = 42 rad/sec, and an emitter T4 at frequency 5.5 MHz. Dashed line is an extrapolation based 
on the flow model. 

FIG. 19. (a) Modulus of the Fourier transform of the vorticity obtained from the data in Fig. 18. 
(b) Normalized spatial correlation function of the vorticity. Dashed line-sound scattering results, 
solid line- calculations based on PIV velocity profile measurements as fig. 5 (with extrapolation up 
to the cell wall) with ro = 14mm and = 42 rad/sec. 

FIG. 20. Temporal dependence of the vortex location detected by PIV (dots in the upper set) 
and by sound scattering (squares in the lower set) measurements. 

FIG. 21. Auto-correlation functions for two sets of the data in Fig.20 as a function of the time 
delay. Solid line- for sound scattering, dash-dotted line-for PIV. 

FIG. 22. (a) Temporal dependence of the peak of the modulus of the velocity structure function 
for the same measurements as in Fig.20. (b) Auto-correlation function for the data in (a) shows 
that the vortex precession frequency is half of the peak frequency. 
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FIG. 23. Vortex precession frequency obtained from the auto-correlation functions of the 
far-field scattering peak time dependence (similar to the plot in Fig.22b), as a function of the 
frequency of the plate rotation. The proportionality factor is 0.22 zt 0.02. 
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